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Abstract. - We study vortex lines in high-temperature superconductors with columnar defects 
produced by heavy ion irradiation. We reconsider scaling theory for the Bose glass transition 
with tilted magnetic fields, and propose, e.g., a new scaling form for the shape of the Bose glass 
phase boundary, which is relevant for experiments. We also consider Monte Carlo simulations 
for a vortex model with a screened interaction. Critical exponents are determined from scaling 
analysis of Monte Carlo data for current- voltage characteristics and other quantities. The 
dynamic critical exponent is found to be z = 4.6 ± 0.3. 



A rich variety of new vortex phase transitions has been established in high-temperature 
superconductors H. Bose glass physics was originally suggested, by Nelson and Vinokur ||, 
to apply for vortex phase transitions in systems with artificially introduced columnar defects, 
produced as permanent damage tracks from heavy ion irradiation of the sample. Columnar 
defects act as optimal pinning centers, leading to a considerable increase of critical currents 
and fields, as compared to the unirradiated sample || |^, m. Bose glass theory makes a 
set of distinct predictions, e.g., for the dynamic scaling of vortex transport properties, and 
the response to tilted magnetic fields 0, H, M. Such behavior has been observed in various 
experiments f| |, |, [[^ @ . 

Vortex lines in 3D superconductors with columnar defects can be mapped to imaginary time 
world lines of bosons in 2+1 dimensions on a disordered substrate The superconducting 
glass phase for the vortex lines corresponds to the insulating Bose glass (BG) phase for the 
bosons. The dissipative vortex line liquid phase corresponds to the superconducting phase 
for bosons. This mapping gives useful information about the equilibrium properties in the 
columnar defect problem p3|. However, dynamical properties do not follow from the mapping. 
A tilted magnetic field H±_ with respect to the columns enters like an imaginary vector potential 
for the bosons leading to a localization problem in non-Hermitian quantum mechanics, 
which has received considerable attention recently |l4|, |l5) . The vortices want to stay localized 
on the columns, which results in a transverse Meissner effect with a divergent tilt modulus C44. 
The BG phase boundary in the (T, ffjj-plane has a sharp cusp at H± = 0, which distinguishes 
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the Bose glass from the isotropic point-defect vortex glass which has a smooth phase boundary 
[pi . Such a feature in the transition line is seen in experiments [^[ |l^| , and provides support 
for the Bose glass theory. Other types of correlated disorder, like splayed columnar defects Jl6| ], 
and planar defects 0, also strongly influence the properties of the glass phase. 

In this paper we reconsider scaling theory for the Bose glass transition for the case of 
tilted magnetic fields, which leads to certain important modifications of predictions in earlier 
work e.g., for the form of the BG phase boundary. We further present Monte Carlo (MC) 
simulations for various thermodynamic and dynamic quantities, which verifies the scaling 
predictions derived below, and give an estimate of the dynamical critical exponent z. 

We first discuss scaling theory. A detailed scaling theory has been developed for the 
superconducting phase transition by Fisher et al. ]l7t , and generalized to anisotropic Bose glass 
scaling in Refs. |^, ^) . In the present problem it is crucial to correctly distinguish between the 
scaling of the B and H fields. These relations have sometimes been mixed up in the literature. 
On approaching the transition the correlation lengths in the directions perpendicular and 
parallel to the columnar defects are assumed to diverge as £ = £x ~ \T — Tbg\ u and £m ~ 
respectively, where £ is an anisotropy exponent. Due to screening of the interaction and 
the correlated disorder, the correlation volume is anisotropic with £ = 2 In case of 

unscreened long-range interactions we expect instead £ rs 1. The correlation time diverges 
as t ~ £ z where z is the dynamical critical exponent. The vector potential enters in the 
combination V — (27ri/<E>o)A, and therefore scales as A±_ ~ A\\ ~ £7 . From hyperscaling 
the free energy density scales as / ~ where d is the dimension, and therefore the 

current density scales as J± — -gf^ ~ £ 2 ~ rf ~'',^|| = -§x^ ~ £, ld , and the electric field as 

E ± = ~ £," (1+z \ E\\ = ~ £~ (C+z) - The scaling of the magnetic field is obtained 

from B = -4tt|£ (or from V x H = ^J), which gives H± ~ <f H\\ ~ f- d ~i , and the flux 
density scales as B± — (V x A)j_ ~ £~ 1_ S-B|i = (V x A)|| ~ £ -2 . The appropriate scaling 
combinations involving the magnetic field are therefore H±^ d ~ 2 and H\\£ d - 3+c , which differs 
from those found, e.g., in Ref. Q). This has several experimentally relevant consequences as 
we will discuss below. The linear resistivity scales as 

pi = z dH - 3 - z pi(H ± z d - 2 ) (i) 

P\\ = f-^-'phHxt*- 2 ), (2) 

and the I-V characteristic as 

E ± = r {1+z) Ei(J^ d+<: - 2 ,H^ d - 2 ) (3) 
E\\ = r iC+z) El(J^ d -\H ± C d - 2 ), (4) 

where p and E are scaling functions. For the magnetic permeability we obtain 

fix = £ d - 3 - C fi±(H ± Z d - 2 ) (5) 
M|| = Z d - 5+C H\\(H^ d - 2 ). (6) 

As an independent check of the correctness of this scaling law we can again use the mapping 
to dirty bosons in (2+l)D. The permeability /j,± corresponds to the superfluid density p s for 
the bosons, whose scaling was derived in Ref. [|l3| with the result p s ~ £ 3_d ~C. This scaling 
law coincides with Eq. for d = 3 Equation (|^) corresponds to the compressibility of 
the bosons k ~ £^ +1 , again with agreement for d = 3. To obtain the scaling of the BG phase 
boundary for small tilt, Eq. (|) gives n ± = t- u{ - d - 3 -^p,±{H 1 t-^ d -^), where t = \T - T BG |, 
and following Ref. we assume that the Bose glass phase persists up to a finite tilt, so p±(x) 
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has a singularity at a finite value x — x c . At this point the critical tilt field vanishes when 
T -> T BG as 



for d = 3. Equations differ from the corresponding ones in Ref. j2|, and in particular our 
Eq. (0) does not include a factor 3 in the exponent. The shape of the phase boundary still has 
a cusp at H± = 0, but not as sharp as the earlier predictions suggested. Another consequence 
of experimental interest is the scaling relation p± ~ for the resistivity at T = Tbg- 

We now turn to our Monte Carlo simulations. We concentrate on low magnetic fields, 
B < <E>o/A 2 , where the vortex interaction is effectively short ranged, and also restrict our 
study to fields below the matching field, where the density of columnar defects is greater 
than the vortex density. Vortex dynamics has been studied previously in simulations of the 
nonlinear current-voltage (I-V) characteristics j^] for finite applied supercurrents. In this paper 
we present extensive simulation results for the linear resistance, and compare with dynamical 
scaling of the I-V characteristics. We choose to study the dirty Boson action with an onsite 
repulsion , which gives a coarse-grained representation of the physical system, 



where the integer variable m = (rn x ,m y ,m z ) is the vorticity vector on the links of a simple 
cubic lattice of size fl = L x L x Lu with lattice constant a = 1. The partition function is 
Z = Tr e~@ H , where the trace is the sum over all vortex line configurations with no open ends, 
thus satisfying the constraint V • m = 0, and p = 1/T is the inverse temperature. Periodic 
boundary conditions are used in all directions to eliminate surface effects. The columnar 
defects are aligned in the ^-direction, and are modeled as a uniformly distributed random site 
energy g(r±), which is constant in the z-direction. The applied magnetic field is included as 
a finite fixed net density of vortex lines. We also consider tilted magnetic fields by allowing 
global vortex line fluctuations in the direction perpendicular to the columns, with a bias term 
— ^ r H x m rx added to the Hamiltonian. Although not realistic in detail, this model should 
fall into the universality class of strongly screened vortices with columnar defects and therefore 
give the critical exponents correctly for low enough magnetic fields. 

Our Monte Carlo method for vortex line models has been described elsewhere and is only 
outlined here jl9[ . The starting configuration is taken as an assembly of straight vortex lines 
penetrating the system in the z-direction. The MC moves are attempts to add closed vortex 
loops around randomly selected elementary plaquettes of the lattice. The trial moves are 
accepted with probability 1/(1 + e l3AH ), where AH is the total change in energy. 

I-V characteristics can be calculated by the following method J|, ^9|. Each time a loop is 
formed it generates a voltage pulse AQ = ±1 perpendicular to its plane, the sign depending on 
the orientation of the loop. This leads to a net electric field E(t) = -^J v (t), (in the following 
we set h/(2e) = 1) where the vortex current density at MC time t is given by J v (t) = ^j, and 
At = 1 for one full sweep through the system, where, on average, an attempt is made to create 
or destroy one vortex loop on every plaquette of the lattice. The nonlinear I-V characteristic 
can be modeled as the electric field E, due to vortex current response in the presence of a 
uniform Lorentz force on the vortex lines, proportional to the applied current density J. The 
linear resistance can be calculated from the equilibrium voltage fluctuations via the Kubo 
formula Q, R = ± J2Z-oo At [(V(*)^( ))]- Here (' ' ") denotes thermal average and [• • ■] 
disorder average. 

To study the response to tilt, we consider an ensemble with fluctuating winding number. 
We achieve this by including global MC moves where vortex lines are inserted across the whole 
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Fig. 1. - MC results for the linear resistivity vs. system size L at the Bose glass transition in the 
directions parallel and perpendicular to the columnar defects. From the slope of the curves the 
dynamical critical exponent is estimated to be z ~ 4.6 for both directions. 



Fig. 2. - Nonlinear current-voltage characteristics at the Bose glass critical point. Dashed lines have 
slope one (Ohmic response). Solid lines have slope corresponding to the value z = 4.6 from the linear 
resistivity in Fig. 1. 



system in the direction perpendicular to the columns. The uniform part of the flux density 
is proportional to the winding number, B v = ^jf- = jjX}r TOr "- ^ n a tilted magnetic field 
[(-B_i_)] is non-zero above the BG temperature and zero below. The magnetic permeability in 
the direction perpendicular to the columns is given by 



_ a [(£„)] _ n 

dH v ~T 



{Bl)-{B u y , (9) 



and is related to the tilt modulus by C44 ~ I/MX- I n the analogy to dirty bosons, fi± 
corresponds to the renormalized superfiuid density p s of the bosons. 

The details of our MC calculations are as follows. We use a net density of vortex lines 
corresponding to filling of / = 1/2 flux quanta per plaquette. We use a disorder potential g 
uniformly distributed in [0, 1]. The model, Eq. (||), has been extensively studied in the context 
of dirty bosons, and has a glass transition at Tbg ~ 0.248 |l3| (in case of no tilted magnetic 
field). In finite systems the scaling functions obtain extra arguments £/L and £n/Ln. To 
enable finite size scaling of MC data we use system sizes L z = cL 1 *, making the ratio Lii/L^ 
a constant, which eliminates one argument from the scaling functions. We use £ = 2, c = 1, 
and lattice sizes L = 4, . . . , 16. For each realization of the columnar disorder, up to 2 • 10 4 MC 
sweeps were discarded for warmup, and measurements taken during up to 2 • 10 6 sweeps. Up 
to 10 3 disorders were used to obtain small statistical errors in the disorder averages. 

We now turn to our MC results, starting with the magnetic field applied along the columns. 
Figure [l] shows a log-log plot of the linear resistivities p± in the direction perpendicular to 
the columnar defects, and pu parallel to the columns, vs. system size L at T = Tbg- Power 
law fits to the data (solid lines) allow z to be determined according to Eqs. (|l],||); which gives 
z — 4.6 ± 0.3 in both directions. 

Figure || shows a finite size scaling data collapse according to Eqs. (j3|^) of the nonlinear 
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Fig. 3. - Finite size scaling of MC data for B± for various sizes L and applied perpendicular fields 
H± at T = Tbg- The dashed line indicates the form B± ~ H^ c , with £ = 2, for H±L > 1. 



Fig. 4. - MC data for magnetic permeability in a tilted magnetic field. The critical temperature 
T C (H±) is roughly estimated as the point where the curves for different L intersect (filled circles in 
the inset). Inset: The resulting BG phase boundary. The solid straight line corresponds to a power 
law T BG - T ~ l^xl 17 " with v = 1.0. 



current-voltage characteristics for currents applied perpendicular and parallel to the columnar 
defects, at T — Tbg- In the limit of small currents we observe Ohmic response when the 
nonlinear response length scale exceeds the system size. In a finite range of sufficiently small 
currents the I-V characteristics show powerdaw behavior as indicated by the straight lines in 
the plot. The solid straight lines are given by the scaling forms in Eqs. @ and @), using the 
value z — 4.6 determined in Fig. [l] from the linear resistivity. For higher currents there are 
clearly visible deviations, and a crossover to different power laws takes place. If this range of 
large currents is fitted to Eqs. one has to assume different values of z in the _L and || 

directions: z± « 6, zn ~ 4, in agreement with Ref. [||. 

We next consider tilted magnetic fields. As before a constant field was applied parallel to 
the columns, corresponding to a half filling, and in addition we apply a small perpendicular 
field H±. We first consider the effect of H± on B± at the Bose glass temperature T = Tqq, 
shown in Fig. [| The data collapse for different system sizes and tilts verifies the scaling 
form B± = L^ 1 ~^B±{H±L d ^ 2 ). The dashed line corresponds to the power law scaling 
form B± ~ H]^~ for scaling variable H±L > 1. To locate the BG phase boundary as 
the magnetic field is tilted, we use Eq. (^|), which, for finite systems and d = 3, gives 
fj,X = L~^jix([T — T C (H±)]L 1 '" ). Here v" is the correlation length exponent for finite tilt 
fields, which belongs to a separate universality class M. At the critical temperature in a finite 
system, this should scale as fix ~ L _< * for large enough L. In Fig. || we use the crossing point 
for MC data for small system sizes as a crude estimate of the critical temperature. However, 
rather large corrections to scaling are expected due to the closeness of the BG fixed point, 
and we were therefore not able to determine v" . The inset in Fig. |] shows the transition 
temperatures (filled circles) vs. H±. The solid straight line represents Eq. (F?j), using the value 
v = 1.0 for the BG transition jl3). Thus the form of the BG phase boundary given by Eq. (Q) 
seems consistent with the simulation. 
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Finally we will compare our results with experiments. Transport measurements for 
Tl2Ba2CaCu208 (Tl-2212) thin films with columnar defects and zero tilt Jll| obtained 
z k 4.9, v k 1.1, C ps 1.9, which closely agree with our exponents. Experiments for Tl-2212 
with tilted magnetic fields |7j] obtained z" « 4.4, v" ~ 1.8, where the double prime indicates 
exponents for the tilted fixed point. A heavily twinned YBCO single crystal, without arti- 
ficial columnar defects, gave fl2fl v{z - 2) w 2.8, for zero tilt, which is close to our number 
v{z — 2) « 2.6. Quite different experimental results were found in Ref. Experiments on 
YBa2Cu3C>7 (YBCO) single crystals with columnar defects in low magnetic fields (0-6.3 kOe) 
H suggest z w 2.2,1/ rs 1,£ rj 1. The isotropic scaling indicates that the vortex interaction 
is not effectively screened p3l, and a model with a longrange vortex interaction may be 
suitable. We have also done some (limited) simulations for a model with unscreened longrange 
interactions, and find a dynamical exponent z f» 2, which is considerably smaller than the one 
for strong screening. Experiments with tilted magnetic fields (lj, ^|, |§| show a sharp cusp in 
the BG phase boundary, in agreement with the Bose glass theory. We note, however, that the 
actual shape of the BG phase boundary in the figures in Refs. jl^, [)| |^], appears to reasonably 
well fit our 7bg(0) - Tbg(H±) ~ \Hj_\ 1/,J for small H±, with exponent l/v « 1.0 instead of 
1/3^ from Ref. g). More experimental data for the precise shape of the BG phase boundary 
and for other quantities would be useful to test our new scaling relations. 

In summary, we have analyzed the Bose glass transition in superconductors with columnar 
defects by means of scaling theory and Monte Carlo simulations. For magnetic fields tilted 
away from the direction of the columns we suggest a form of the BG phase boundary for 
small tilts, which appears consistent with our simulations and with experiments. Our value 
z = 4.6 ± 0.3 for the dynamical critical exponent at the BG transition in a strongly screened 
model is in good agreement with recent experiments. 

We gratefully acknowledge very useful discussions with D. R. Nelson, S. M. Girvin, and S. 
Teitel. This work was supported by the Swedish Natural Science Research Council, by the 
Swedish Foundation for Strategic Research (SSF), and by Parallelldatorcentrum (PDC), Royal 
Institute of Technology. 
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